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ABSTRACT 

Complex organic molecules have previously been discovered in solar type protostars, raising the 
questions of where and how they form in the envelope. Possible formation mechanisms include grain 
mantle evaporation, interaction of the outflow with its surroundings or the impact of UV/X-rays inside 
the cavities. In this Letter we present the first interferometric observations of two complex molecules, 
CH 3 CN and HCOOCH3, towards the solar type protostar IRAS16293-2422. The images show that 
the emission originates from two compact regions centered on the two components of the binary 
system. We discuss how these results favor the grain mantle evaporation scenario and we investigate 
the implications of these observations for the chemical composition and physical and dynamical state 
of the two components. 

Subject headings: ISM: abundances — ISM: individual (IRAS16293) — ISM: molecules — stars: 
formation 



1. INTRODUCTION 

Solar type Class protostars are characterized by be- 
ing observable only at millimeter to far infrared wave- 
lengths. This is because the central forming stars are 
surrounded and obscured by massive envelopes, making 
them the coldest known protostars. Yet, the cold en- 
velopes hide inner warm regions where the grain mantles, 
built during the pre-collapse phase, evaporate. Follow- 
ing gas phase reactions between the evaporated mantle 
molecules (also called first generation molecules) , organic 
complex molecules (second generation molecules) form in 
these regions. While this picture has now been widely ac- 
cepted for massive protostars (see e.g. Kurtz et al. 2000; 
van Dishoeck & Blake 1998), it is only recently that evi- 
dences in support of parts of this scenario have been accu- 
mulated for low mass protostars. Actually, so far the ex- 
istence of inner warm regions with evaporated mantles in 
solar type protostars has been argued indirectly from the 
sophisticated analysis of molecular multi-frequency sin- 
gle dish observations (Ceccarelli et al. 2000a, b; Schoier 
et al. 2002; Maret et al. 2002, 2004; Doty et al. 2004). 
A survey on almost a dozen of Class sources has shown 
that, typically, the inner warm regions have sizes of a few 
tens of AUs (Maret et al. 2004), i.e the size of the Solar 
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Nebula. Now that the existence of these regions has been 
demonstrated, the questions of their nature and chemi- 
cal composition have become relevant because of the link 
with the formation of our own Solar System. What is the 
inventory of organic molecules of such regions? What 
are their origin and evolution? In particular, could these 
molecules be incorporated into the planet-forming disks 
surrounding the protostars? 

Based on theoretical arguments, no complex, second 
generation molecules should have the time to form in so- 
lar type protostars. The reason is that the gas crossing 
time of the warm regions is estimated to be much shorter 

10 3 y r )> than the predicted chemical formation time 

10 4 yr; Schoier et al. 2002). But our theories have 
somewhere flaws, since complex organic molecules have 
been detected in the two (out of two) sources where they 
have been searched for: IRAS16293-2422 (Cazaux et al. 
2003) and NGC1333-IRAS4A (Bottinelli et al. 2004). 
We therefore may even doubt the basic prediction that 
those molecules originate in the warm inner regions of the 
envelopes! It is, for example, possible that the process 
which is responsible for the release in the gas phase of 
the mantle constituents, is not the thermal evaporation, 
but the interaction of the outflow with the surround- 
ings, or even of UV/X-rays with the cavities excavated 
by the outflow (Schoier et al. 2002, 2004). Alterna- 
tively, the emission could come from disk surfaces which 
would have been heated by accretion shocks (Schoier et 
al. 2004). The proof that the complex organic molecules 
detected by Cazaux et al. (2003) and Bottinelli et al. 
(2004) originate from the inner warm regions (envelope 
or disk), called hot corinos (Ceccarelli 2004; Bottinelli et 
al. 2004), was still missing until today. 

In this Letter we bring the missing proof and report 
interferometric observations of two complex, second gen- 
eration molecules, methyl cyanide and methyl formate, 
towards the prototype of the hot corinos, IRAS16293- 
2422 (hereafter IRAS16293). 

2. OBSERVATIONS AND RESULTS 
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Observations of IRAS16293 (a(2000) = 16 h 32 m 22!6, 
(5(2000) = -24°28'33") were carried out at the IRAM 
Plateau de Bure Interferometer on February 1 st and 
March 25 th 2004 in the B and C configurations of the 
array. Five CH3CN transitions at 110.4 GHz and 4 
HCOOCH3 transitions at 226.4 GHz were obtained si- 
multaneously, along with the continuum emission at 
3 mm and 1.3 mm. The receivers were tuned single 
side band at 3 mm and double side band at 1.3 mm. 
CH3CN and HCOOCH3 transitions were covered with 
two correlator units, each of 40 and 80 MHz bandwidth 
respectively. Typical system temperature were 250 K 
(USB) at 3 mm and 500 K (DSB) at 1.3 mm. Phase and 
amplitude calibrations were obtained by observing the 
nearby point sources 1514-241 and NRAO 530 every 20 
minutes. Bandpass calibration was carried out on 3C273 
and 0851+202 and the absolute flux density scale was 
derived from MWC349, 3C345 and 0923+392. Data cal- 
ibration was performed in the antenna based manner and 
uncertainties are less than 10% at 3 mm and less than 
20% at 1.3 mm. Flux densities were obtained from vis- 
ibilities using standard IRAM procedures. Continuum 
images were produced by averaging line-free channels. 
Line maps were obtained by cleaning line images after 
subtraction of the continuum directly from the visibili- 
ties. 

Figures ^a and Q}c show the integrated line emission 
of CH 3 CN and HCOOCH3, averaged over all the tran- 
sitions listed in Table 1 for each molecule. Note that 
the energy of the upper level of the CH3CN transitions 
decreases with frequency, but that maps averaged over 
each individual transition do not show any significant 
difference. This means that the emitting region does not 
depend on the energy of the transition, i.e. on the exci- 
tation conditions (but rather on a jump of the molecular 
abundances). Continuum emission at 3 and 1.3 mm is 
displayed in Figures ^-b and ^d respectively . These 
maps show two components which are spatially coinci- 
dent with the centimeter wavelength emission regions A 
and B mapped by Wooten (1989) and with the millime- 
ter wavelength emission regions MM1 and MM2 mapped 
by Mundy et al. (1990, 1992). As already noted in the 
previously mentioned works (see also Looney, Mundy & 
Welch 2000; Schoier et al. 2004), the south-east region 
("source A" or "MM1") is the weakest in the contin- 
uum but brightest in line emission. On the contrary, the 
north-west region ("source B" or "MM2") is the bright- 
est in the continuum and weakest in line emission. 

Table 1 gives the intensities and sizes of the line and 
continuum emissions. Within the errors, we recover all 
the line emission measured by the IRAM 30m (Cazaux 
et al. 2003). 

Finally, the spectra at 1.3 mm of emission regions A 
and B are shown in FigureEl assuming a Vlsr=3.9 (Fig- 
ure 0-a) and 2.7 km s _1 (0-b) respectively (see discus- 
sion in the next section). Again, note that the sum of A 
and B (Figure 0-c) reproduces very well the features in 
the spectrum obtained at the IRAM 30m (Cazaux et al. 
2003). 

3. DISCUSSION AND CONCLUSIONS 

The most important result of the presented observa- 
tions is that the complex molecules observed by Cazaux 
et al. (2003) originate in two compact regions, whose di- 



ameters are about 1.5" (A) or less (B), as shown in Fig- 
ure H (and Table 1). This goes along with the fact that 
the images do not show any evidence of emission asso- 
ciated with the molecular outflows seen at larger scales 
(see below for the discussion on the line profiles). Also, 
following the remark in the Introduction, the Plateau de 
Bure images do not reveal any evidence of emission from 
cavities excavated from the outflows either. Indeed, the 
cavities are typically ~ 20" in size, located ~ 30" away 
from the source along the 12 CO outflow (see e.g. Arce & 
Sargent 2004), which would have been easily detected by 
the Plateau de Bure. The two regions where the molecu- 
lar emission comes from, are compact, and while source A 
is barely resolved in the 1mm images 8 , source B is unre- 
solved. The measured sizes (1.5" at 160 pc correspond to 
a radius of about 120 AU) are remarkably consistent with 
the emission coming from a region where the dust tem- 
perature exceeds 100 K in source A (150-200 AU: based 
on multi-frequency single dish observations: Ceccarelli et 
al. 2000a, b; Shoier et al. 2002), and therefore, where the 
grain mantles evaporate. Thus, these observations sup- 
port the basic prediction (from the modeling of the single 
dish observations; Ceccarelli et al. 2000a,b; Cazaux et al. 
2003) that a hot corino with a radius of about 150 AU 
exists inside the cold envelope of IRAS16293, and that 
in that region, complex molecules are formed because of 
grain mantle evaporation. 

In addition to that, the Plateau de Bure observations 
confirm that the two sources A and B are different, as 
noted by previous authors (Wootten 1989; Mundy et al. 
1990, 1992). They differ in line intensities and extent 
(Fig. and this corresponds to a difference in their 
chemical composition. But, before discussing this point, 
it is necessary to address the second most striking dif- 
ference in the two sources: their line profiles (Fig. 
Source A has clearly broadened spectra (FWHM ~ 8 
km s _1 ), while source B shows apparently much narrower 
profiles (FWHM ~ 2 km s _1 ). Furthermore, the lines of 
source B seem to peak at Vlsr=2-7 km s" 1 , whereas the 
parent cloud velocity is at Vlsr=3-9 km s^ 1 (compatible 
with the spectra of source A, although given the broad 
profiles it is difficult to precisely determine the Vlsr of 
source A). Note that the cloud's Vlsr (3-9 km s _1 ) is 
very nicely measured by the two CN absorption lines in 
the observed band 9 (Figure As said, no evidence 
of outflowing gas is seen in the images, and also the 
broad line profiles of source A are consistent with gas in- 
falling towards a ~1 M Q object (Ceccarelli et al. 2000a, b; 
Schoier et al. 2002). Therefore, both the images and 
line profiles of source A are fully consistent with the hot 
corino hypothesis. 

The case of source B is less obvious: why does this 
source have narrower lines and why do they peak at 2.7 
km s _1 ? But is this true? A very careful look at the B 
spectra rises doubts. Indeed, all the B lines have a sec- 
ond small peak — sometimes at the limit of the noise — 
on the red-shifted side of the spectrum. This second 

8 It is likely that the emission from source A originates in the disk 
around this source. A detailed study of this aspect is postponed 
to a forthcoming paper. 

9 The absorption originates in the foreground (envelope + cloud) 
cold gas, which absorbs the photons emitted in the hot corino re- 
gions (the CN emission component, being extended, is filtered out 
by the interferometric observations). 
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Fig. 1. — (a) Line map of CH3CN, averaged over the five transitions listed in Table 1; the rms is 3 mjy beam -1 and contours range from 
15 to 150 mjy beam -1 in steps of 15 mjy beam -1 , (b) Continuum emission at 3 mm, with an rms of 3 mjy beam -1 ; contour levels range 
from 20 to 220 mjy beam -1 in steps of 20 mjy beam -1 , (c) Line map of HCOOCH3-A and -E, averaged over the four transitions listed 
in Table 1; the rms is 8 mjy beam -1 and contour levels range from 20 to 200 mjy beam -1 in steps of 20 mjy beam -1 , (d) Continuum 
emission at 1.3 mm, with an rms of 15 mjy beam -1 ; contour levels range from 50 to 600 mjy beam -1 in steps of 50 mjy beam -1 . Beam 
sizes are 4'.'7 X l'/6 and 2'.'2 X 0'.'9 at 3 and 1.3 mm respectively. 



TABLE 1 
Line and continuum emission 



Molecule 


Transition 


E up a 


Frequency (GHz) 


(Integrated) intensity 


Size (") b 




ARA.ADec (") c 








(cm -1 ) 


or Wavelength 


A 


B 


A 


B 


A 


B 


CH 3 CN 


65,0 


- 55,o 


137.1 


110.330 


1.39(0.09) 


0.51(0.06) 


0.8(0.2) 


< 0.8 


3.6,-3.3 


-0.0,1.1 




64,0 


— 54,0 


92.4 


110.350 


1.70(0.08) 


0.38(0.06) 


0.8(0.2) 


< 0.8 


3.7,-3.2 


-0.1,0.7 




63,0 


- 53,0 


57.6 


110.364 


2.75(0.08) 


0.67(0.05) 


0.9(0.1) 


< 0.8 


3.7,-3.2 


+0.1,0.6 




62,0 


— 52,0 


32.8 


110.375 


2.56(0.08) 


0.51(0.05) 


1.2(0.1) 


< 0.8 


3.6,-3.3 


+0.1,0.9 




61/0,0 


- 5i/o,o 


17.9/12.9 


110.381/110.383 


5.26(0.10) 


1.05(0.07) 


0.9(0.1) 


< 0.8 


3.6,-3.3 


-0.0,0.6 


HCOOCH3-E/A 


20 2 ,i9 


- 192,18 


83.5 


226.713/226.718 


6.44(0.41) 


3.36(0.36) 


1. 4(0. 2)x0. 7(0.1) 


0.8(0.2) 


3.6,-3.5 


-0.1,0.7 




20i,i9 


- 19l,18 


83.5 


226.773/226.778 


11.13(0.49) 


4.05(0.39) 


1. 6(0. 1)x0. 8(0.1) 


0.8(0.1) 


3.7,-3.3 


-0.0,0.7 


continuum 








3 mm 


0.17 


0.26 


3.4x1.4 


1.5x0.8 


3.5,-2.9 


+0.2,0.7 


continuum 








1.3 mm 


0.77 


1.02 


3.7x1.2 


1.5x0.8 


3.5,-3.3 


+0.2,0.6 



Note. — The line intensity is in units of Jy km s 1 , while the continuum is in Jy. Errors are given in parentheses. 



a Encrgy of the upper level of the transition. 

b FWHM of gaussian fit. 

c Rclativc peak position from the fit. 
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1 *0 f (a) PdB,' F LM '= 2.7 km/s 




226.79 .77 .76 .72 .71 .70 .67 .66 .65 
Rest Frequency (GHz) 

Fig. 2. — (a) Spectra at 1 mm averaged over the emission re- 
gions of sources A (dotted line) and B (solid line), and displayed 
for Vlsr=2.7 km s — 1 . The velocity resolution is 0.4 km s" 1 and 
rms are 0.08 and 0.06 Jy for A and B respectively, (b) Same 
as (a) but with Vlsr = 3-9 km s — 1 . (c) The solid line represents 
the spectrum obtained at the IRAM 30m (Cazaux et al. 2003) 
with a velocity resolution of 1.3 km s and rms of 0.11 Jy. Over- 
layed in dotted line is the sum of the Plateau dc Bure spectra of 
sources A and B (rms=0.10 Jy). Vertical lines indicate the fre- 
quencies of HCOOCH3 (A and E) lines (emission) and CN lines 
(absorption). Right to left are CN 235 — li 3, 231 — li 1, 

2 3 3 - 1 1 1 , HCOOCH3-E and -A 2O2 19 - 19 2 is, HCOOCH3-E 

2 ' 2 2 1 2 
and -A 20i,i9 - 19l.l8- 



peak could indeed be part of the line itself, which would 
be strongly self-absorbed at Vlsr=3.9 km s -1 . If this 
is the case, the linewidths of source B would be ~ 4-6 
km s" 1 (Fig. EJl, similar to the linewidths measured to- 
wards source A. Note that the blue peak is expected to be 
brighter than the red peak in the case of optically thick 
lines from infalling gas (Leung & Brown 1977; Zhou 1992, 
1995; Choi et al. 1995), so it would be consistent with 
the ~ 1 Mq hot corino hypothesis of source B too. This 
alternative explanation, optically thick lines in source B, 
is therefore very appealing and worth exploring in some 
detail. Using the LTE approximation, the required col- 
umn densities for the CH3CN and HCOOCH3 lines to be 
optically thick are N ~ 10 16 cm -2 and N ~ 10 17 cm" 2 
respectively. These values are about one order of mag- 
nitude larger than the column density derived in B from 
the emission lines, assuming that the lines are optically 
thin, LTE populated and that the emission region fills 
up the Plateau de Bure synthesized beam (see below). 
Considering that the three adopted assumptions all un- 
derestimate the true column density, it is indeed possible, 
but not firmly established, that the lines in source B are 
optically thick. Unfortunately, "physical" considerations 



do not help either to distinguish between the two possi- 
ble interpretations, optically thin or thick lines in source 
B. In the first case (B has VLsr=2.7 km s" 1 and FWHM 
~ 2 km s" 1 ), source B would be less massive than A and 
would revolve around it at 1.2 km s" 1 (multiplied by the 
inclination of the orbit), at a distance of 800 AU, which 
is fully consistent with Ma ~ 1-M© (unless the orbit is in 
the sky plane). In the second case (B also has Vlsr=3.9 
km s" 1 and FWHM ~ 4-6 km s _1 ), A and B have com- 
parable masses (similar FWHM), but B is more compact. 
Future high resolution observations of optically thin lines 
are hence required to definitely settle the question. 

As said, the nature of source B affects the determi- 
nation of the molecular abundances in this source, and 
hence, how much the chemical composition of the A and 
B hot corinos differ. Using the relation between column 
density and observed continuum flux density, and using a 
dust opacity k v of 0.8 cm 2 g _1 at 1.3 mm and a dust tem- 
perature Td of 40 K (e.g. Walker et al. 1986), we derive 
molecular hydrogen column densities from the 1.3 mm 
continuum emission equal to N(H.2, A) = 3.5 x 10 24 cm -2 
in A and A(H 2 ,B) = 1.7 x 10 25 cm" 2 in B (consistent 
with the values of Mundy et al. (1992) 10 ). Using these 
values and the CH3CN total column densities derived 
from the rotational diagram method (N\ = 1-7 x 10 15 
cm" 2 and N B = 3.9 x 10 14 cm" 2 ), we get CH 3 CN 
abundances of 4.8 x 10" 10 and 2.3 x 10" 11 for A and 
B respectively, i.e. CH3CN is ~ 20 times more abun- 
dant in A than in B. A rotational diagram could not 
be drawn for the HCOOCH3 transitions as they have 
similar upper energy levels, but assuming T rot ^60 K 
(Cazaux et al. 2003) we get HCOOCH3 column densi- 
ties of N A = 2.6 x 10 16 cm" 2 and N B = 8.4 x 10 15 cm" 2 , 
i.e HCOOCH3 abundances of 7.5 x 10" 9 and 4.9 x 10" 10 
for A and B respectively, a factor 15 difference 11 . Note, 
however, that since B is unresolved in the line emission, 
and the lines could be optically thick (see above discus- 
sion), the molecular abundances quoted for B may be 
underestimated by about an order of magnitude (note 
also that the region of molecular emission could be more 
compact than the continuum emission region). There- 
fore, resolving the problem of the nature of the observed 
line emission in source B is crucial, not only to deter- 
mine the dynamical state of this source, but also to cor- 
rectly assess the difference in the chemical composition 
of sources A and B. 

In summary, these interferometric observations with 
the Plateau de Bure show unambiguously that the ob- 
served complex organic molecules are emitted from two 
compact (< 1.5") regions, which is consistent with grain 
mantle evaporation. No indication was found in support 
of other formation mechanisms for these molecules. Two 
scenarios have been proposed to explain the observed dif- 
ferences in sources A and B. However, the available data 
are insufficient to fully understand how and why the two 
components differ in their dynamical and chemical state, 
and more interferometric observations (e.g. of optically 
thin lines) are needed in order to solve the issues raised 
in this Letter. 

10 We postpone the detailed analysis of the continuum emission 
to a forthcoming paper. 

11 Note that the abundances quoted in Cazaux et al. differ 
from those derived here because of the different estimate of the H2 
column density. 
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